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Abstract 

In this article, we consider a small 50 (2)-equi variant perturbation of a reaction-diffusion 
system on the sphere, which is equivariant with respect to the group 50(3) of all rigid rotations. 
We consider a normally hyperbolic 50 (3)-group orbit of a rotating wave on the sphere that 
persists to a normally hyperbolic 50(2)-invariant manifold M(e). We investigate the effects 
of this forced symmetry breaking by studying the perturbed dynamics induced on M(e) by 
the above reaction-diffusion system. We prove that depending on the frequency vectors of the 
rotating waves that form the relative equilibrium 50(3)ito, these rotating waves will give 50(2)- 
orbits of rotating waves or 50(2)-orbits of modulated rotating waves (if some transversality 
conditions hold). The orbital stability of these solutions is established as well. 

1 Introduction 

In mathematics and physics the phrase "symmetry breaking" has distinct meanings. The first refers 
to the frequently observed phenomenon that a configuration of a physical system satisfying a law 
(a set of equations) which is invariant under a group of transformations, may itself be invariant 
under a subgroup of this group. This is referred to as spontaneous symmetry breaking. The second 
meaning refers to the problem of explicitly adding symmetry breaking terms to the equations which 
describe the system. This is called induced or forced symmetry breaking. In another words, forced 
symmetry breaking means the systematic study of an equivariant system of differential equations 
which is perturbed slightly so that it loses its symmetry properties partially or completely. 
The main motivation of this article is the presence of spiral waves in excitable media, especially in 
cardiac tissue. Spiral waves arise as stable spatio-temporal patterns in various chemical, physical 
systems and biological systems, as well as numerical simulations of reaction-diffusion systems on 
excitable media with various geometries. Excitable media are extended non-equilibrium systems 
having a uniform rest state that is linearly stable but susceptible to finite perturbations. Spiral 
waves have been observed experimentally, for instance, in catalysis of platinum surfaces, Belousov- 
Zhabotinsky chemical reactions ^!5]> Rayleigh-Bernard convection, and the most important, cardiac 
tissue 9 . Numerical simulations of spiral waves have been done, for example in [HI2irn H I35 1 l3fi l l?F7]. 
Winfree |31| found spiral waves in the Belousov-Zhabotinsky reaction. 
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It is now believed that spiral and scroll waves that appear in the heart muscle can lead to cardiac 
arrhythmias (abnormal rhythms in the heart), giving rise to atrial fluttering or ventricular fibril- 
lation. In normal hearts cardiac arrhythmias are rare, but in diseased hearts cardiac arrhythmias 
can become more common. For example, if chambers of the heart become abnormally large, they 
are susceptible to serious arrhythmias in which waves are believed to circulate in a fashion that is 
similar to the circulation of the Belousov-Zhabotinsky waves in a chemical medium. Real human 
hearts are enormously complex three-dimensional structures. In this article, we assume that the 
geometry of the excitable media is a sphere, but the forced symmetry breaking from SO (3) to 
SO (2) allows us to consider either the case of the sphere with a localized inhomogeneity, or the 
case of a domain with SO(2) symmetry, obtained by slightly deforming the sphere. In the case of 
the cardiac tissue, this is clearly an approximation. 

In the planar case, a rigidly rotating spiral wave is an example of wave pattern rotating around a 
center and being well approximated by an Archimedean spiral wave far from rotation center. Near 
the rotation center, there is a core region of the spiral wave, where the front of the wave has a tip, 
whose structure is considered to be the most important in understanding the behavior of the whole 
spiral wave. Typical motions which have been observed for planar spiral waves are rigid rotation, 
quasi-periodic-meandering and linear meandering. Using an ad hoc model, Barkley [2] was the 
first to realize the key importance of the group SE(2) of all planar translations and rotations in 
describing the dynamics and bifurcations of planar spiral waves. From a mathematical point of 
view, for the planar case, rigidly rotating spiral waves are examples of rotating waves, meandering 
spiral waves are examples of modulated rotating waves and linearly drifting spiral waves are ex- 
amples of modulated travelling waves. The first rigorous mathematical theory of the planar spiral 
waves was done by Wulff 32 using Liapunov- Schmidt reduction on scales of Banach spaces. These 
results were generalized by Wulff to general non-compact Lie groups in [SB]- Comanici [0] and 
Chan 4 have independently expanded this work to the rotation group G = SO(3), using different 
approaches. Sandstede, Scheel and Wulff [2Hj proved a finite-dimensional center bundle reduction 
theorem near a relative equilibrium Guq of an infinite-dimensional vector field on a Banach space 
X on which acts a finite-dimensional Lie group (not necessarily compact). Using the results of |28j . 
the Hopf bifurcation from one-armed and multi-armed rotating spiral wave to meandering spiral 
waves can be studied ^2j [2Hj . 

In summary, the Euclidean-equivariant center bundle approach has been remarkably successful at 
explaining many of the experimentally observed dynamics and bifurcations of planar spiral waves. 
However, the Euclidean symmetry is an approximation, since no physical experiment is infinite in 
spatial extent. 

There are indeed experiments on planar spiral waves which exhibit some phenomena which cannot 
be explained by Euclidean symmetry alone: boundary drifting (341 136j . spiral anchoring on localized 
inhomogeneities OEH] and repelling by localized inhomogeneities |25| . Boundary drifting has been 
observed in situations where size of the spiral wave core is comparable to the size of the spatial 
domain of the experiment. In this case, the spiral wave is attracted to the boundary of the domain 
and then drifts around the boundary in a meandering motion. The spiral waves with smaller cores 
do not drift around the boundary: the spiral wave tip goes straight through the boundary without 
any impediment [23] ■ When local inhomogeneities are present, the spiral wave is attracted to these 
inhomogeneities, migrate towards them, and then rotate around them. Some spiral waves are also 
repelled by these inhomogeneities |25| . 

There are also experiments in which two dynamic attractors, called entrainment attractor and 
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resonance attractor are seen, which attractor is observed depends on the initial conditions of the 
experiment JHJ- All of these phenomena break the Euclidean symmetry and as it has been shown 
in LeBlanc and Wulff [2J , these phenomena are generic consequences of imperfect Euclidean sym- 
metry. Specifically, they have studied the effects of translation symmetry breaking on normally 
hyperbolic relative equilibria and normally hyperbolic relative periodic orbits in general systems 
of Euclidean-equivariant differential equations which undergo a small perturbation that breaks the 
translational symmetry, while preserving rotational symmetry. 

Also, both numerical and experimental work suggest that anisotropy also can lead to certain dy- 
namical states for spiral waves which are inconsistent with Euclidean symmetry. Anisotropy can 
lead to phase-locked two- frequency epicycle spiral motions, and to complicated quasi-periodic me- 
andering patterns which have overall discrete rotational symmetries. LeBlanc [221 nas investigated 
the effects the forced rotational symmetry breaking on spiral wave dynamics. 
All previous results are valid for planar spiral waves. The interest to consider spiral waves on non- 
planar surfaces is motivated by the applicability to problems in physiology (cardiology), biology 
and chemistry. Therefore, the study of spiral waves by experiments and numerical simulations of 
reaction-diffusion systems on the sphere and curved surfaces have recently been undertaken. In the 
case of spiral waves on a sphere, the dynamics is expected to be quite different because any spiral 
wave starting from a rotating center cannot end at a point. The number of tips of a wave front 
cannot be odd, and therefore, the dynamics of spiral waves may acquire new features qualitatively 
different from the planar case. 

The dynamics of spiral waves in an excitable reaction-diffusion systems on a sphere was numerically 
investigated by El ES1 GUI > an d [HE] who employ a spectral method using spherical harmonics 
as basis functions. Amjadi treated the case of a periodic oscillating sphere as well. Maselko 
|24| . as well as Maselko and Showalter, performed experiments with Belousov-Zhabotinsky chemical 
waves propagating on the surface of a sphere. The influence of the topological constraints and the 
inhomogeneity in the excitability on the geometry and dynamics of spiral waves on a thin spherical 
shell of excitable media are presented in [H]. Also, rigidly rotating waves on spherical domains have 
been studied using kinematical theory, but we do not intend to describe this here. 
In this article we have studied only the trivial isotropy case, because as far as we are aware, there 
have been observed no m-spiral waves (m > 1) on spherical surfaces. Also, in [Sj it was numerically 
verified that there is a critical size of the sphere below which self-sustained spiral waves cannot 
exist. Therefore, we use a sphere of an arbitrary, but fixed radius r. Guyard and Lauterbach in 
|14| I15| studied forced symmetry breaking perturbations for periodic solutions, but their methods 
applied to the forced symmetry breaking from 5*0(3) to SO (2) do not give the results obtained in 
this article. 

There are almost no experiments or numerical simulations on spherical domains with localized in- 
homogeneities and on approximatively spherical domains, showing phenomena similar to boundary 
drifting, spiral anchoring at localized inhomogeneities or repelling from localized in homogeneities. 
Following the approach for planar spiral waves, this can be investigated and we intend to present 
somewhere else, the numerical simulations illustrating the possibilities of these phenomena for spiral 
waves on spherical domains with localized inhomogeneities and on approximatively spherical do- 
mains. Since for planar spiral waves, forced symmetry breaking from SE(2) to SO (2) was successful 
in explaining phenomena like boundary drifting, spiral anchoring on localized inhomogeneities or 
repelling by localized inhomogeneities, we will treat the same problem here, that is the forced sym- 
metry breaking from SO (3) to SO (2) for rotating waves on the sphere. 
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We consider the group SO (2) as being diffeomorphic with the subgroup of SO (3) defined by 
{e® 9 | 6 G [0, 2tv)}, where Q G so(3) such that \Q\ = 1. Section El is concerned with the set- 
up of the problem of forced symmetry-breaking from SO(3) to SO (2) for a normally hyperbolic 
relative equilibrium SO(3)uq. This is done by studying SO (2)-equi variant reaction-diffusion sys- 
tems on the sphere, that are small perturbations of 50(3)-equivariant reaction-diffusion systems 
on the sphere, having a relative equilibrium SO(3)uo that persists to a normally hyperbolic SO(2)- 
invariant manifold M(e), which is S'0(2)-equivariant diffeomorphic to SO(3). We recall the fact 
that the functional-analytical framework can be found in |33| for the general Lie groups. In Section 
13 we obtain the S'0(2)-equi variant finite-dimensional center manifold reduction and the corre- 
sponding reduced differential equations for the perturbed reaction-diffusion system, which include 
the general form of the S'0(2)-equivariant perturbation on the center manifold near SO(3)uq, 
therefore on M(e). The orbit space reduction methods are presented in Section 0] and we use them 
to project the SO (2)-equivariant perturbed differential equations on M{e) onto the orbit space 
SO (3)/ SO (2), which is diffeomorphic to the unit sphere S 2 . Then, we show that the study of the 
perturbed dynamics (including the orbital stability of the solutions that persist) on M(e) reduces 
to the analysis of some differential equations on the unit sphere S 2 . For e = 0, the dynamics on 
M(e) consist only of rotating waves. Depending on the relation between the frequency vectors of 
these rotating waves and Q , we obtain that the rotating waves will project on the unit sphere S 2 
either onto two antipodal hyperbolic equilibria or onto periodic solutions. Using Implicit Function 
Theorem and the Poincare map, we prove the existence of two smooth branches of hyperbolic equi- 
libria (generically) , and the existence of smooth branches of periodic solutions if some transversality 
conditions hold, and determine the stability of these solutions. Section presents the effects of 
forced symmetry breaking for the rotating waves in the phase space. Namely, depending on the 
frequency vectors of the rotating waves that form the relative equilibrium SO(3)uq, these rotating 
waves will give »SO(2)-orbits of rotating waves or S'0(2)-orbits of modulated rotating waves (if some 
transversality conditions hold). The orbital stability of these solutions is established as well. All 
theorems, for which a reference was not given, are proved in Section |H1 Numerical simulations will 
be presented somewhere else. 



2 Perturbed React ion- Diffusion Systems on the Sphere 

We consider a perturbed reaction-diffusion system of the form 

— (t, x) = DA s u(t, x) + F(u(t, x), A) + eG(u(t, x),x, e, A) on rS 2 , (2.1) 

where r > 0, rS 2 is the sphere or radius r, S 2 is the unit sphere in R 3 and u = (u%, u 2 , 

f d x ... \ 

. . . , un) '■ rS 2 — > M> N for N > 1, As is Laplace-Beltrami operator on rS 2 , D = '■ ' ■ . '■ 

V ... d N j 

with di > for i = 1, 2, . . . , N are the diffusion coefficients, F = (Fx, F 2 , . . . , F N ) : R N x E -> R N 
and G = (G\, G 2 , ... , Gn) ■ ^ N x ?~S 2 x [0, oo) x R —> R N are sufficiently smooth functions for 
|A|, e > small parameters. In the following sections we suppress the parameter A. Therefore, 
let F: R N — > R N and G: R N x rS 2 x [0,oo) -> R^ be sufficiently smooth functions such that 
F(0) = and (7(0, x, e) = for e > small. Using |HH1, the reaction-diffusion system (j2.1|) de- 
fines a sufficiently smooth parameter-dependent local semiflow We recall the notion of linear 
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action of 50(3) on the usual function spaces, that is the function T: 50(3) xY^Y, defined by 
T(A, u)(x) = A • u(x) = u(A~ x x), where A G 50(3), u G Y, x G rS 2 . 

We consider the group 5*0(2) as being diffeomorphic with the subgroup of 50(3) defined by 
{e® e | 6* G [0, 2-71")}, where Q G so(3) such that \Q\ = 1. Thus, we consider that the Lie alge- 
bra so(2) of 50(2) is isomorphic to {Qx \ x G K}. We will study the reaction-diffusion system 
(|2~Tj) on the function space Y = L 2 (rS 2 ,R 7V ). 

Suppose that the superposition operator associated to G, denoted by G: Y a x [0, oo) — > Y, 
0(u, e)(x) = G(u(x),x,e) is 50 (2)-equivariant with respect to the action T restricted to 50(2), 
but not 50 (3)-equivariant with respect to the action T for any e > 0. Then, the local semiflow $ 
is 50 (3)-equivariant for e = 0, but it is only 50 (2)-equivariant for e > small (see J5}). 

Example 2.1. For any x G rS 2 , we denote by r(x) = d(x,l) 2 , where I is the line containing 
the vector Q and d is the Euclidean distance in M 3 . Let us define G(u,x) = G(u,r(x)), where 
G : x R — > is a sufficiently smooth function. Then the superposition operator associated to 
G is SO(2)-equivariant, but is not SO(3)-equivariant in general. 

The definitions of relative equilibrium, relative periodic orbits, rotating waves, modulated rotat- 
ing waves and tip position function for an SO (2)-equivariant or SO (3)-equivariant reaction-diffusion 
system can be found in ||2| for the general Lie groups or in [3] for the rotation group G = 50(3). 
Throughout this paper, we use the following property .26]: AXA^ 1 = B if and only if B = 
AX, where A, B G 50(3) and X G so(3). We also use the following notations: for X = 
a —b \ / c 



-a c I G ao(3), X = I 6 | and |X| 

6 -c 



vV + b 2 + c 2 and IIXII = ^2 



X 



3 The General Form of the 50(2)-Equivariant Perturbation on the 
Center Manifold 

Throughout this paper we suppose that for e = we have a normally hyperbolic relative equilibrium 
uq G Y a such that E uo = 1%. Let <fr(i, «o>0) = e x °'uo- We could have presented only Theorem 13,41 
here, but we have chosen to discuss the SO (2)-equivariant center manifold reduction in Theorem 
21 and Proposition 13.31 for the sake of completion. 



Proposition 3.1 ([51 1111 D-7J). There exists a sufficiently smooth parameter- dependent 50(2)- 
invariant normally hyperbolic manifold M{e) for the reaction- diffusion system \2.1\) such that 
M(0) = 50(3)uo, and, for e > small, M(e) is diffeomorphic to SO(3)uq and this diffeomor- 
phism is SO(2)-equivariant. 

Since S uo = ^3, we have that SO(3)uo is diffeomorphic to 50(3), so M(e) is 50 (2)-equi variant 
diffeomorphic to 50(3). Using UHl HZ1 Ell HH1 we obtain: 

Theorem 3.2. [5] Let L be the linearization of the right-hand side of \2.1\) with respect to the 
rotating wave <fr(i, Uq,0) = e Xot uo in the co-rotating frame, that is 

L = DA S + D u F(u ) - X . 

Suppose that: 
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1. a(L)n{z G C | Re (z) > 0} is a spectral set with spectral projection P* ; and dim(R(P if )) < oo; 

2. the semigroup e Lt satisfies \e Lt \R(i-p„) | < Ce~@ ot for some (3q > and C > 0. 

iei 6e the orthogonal complement of T Uq (SO(3)uq) in E cu = R{P if ). Then, there exists a 
sufficiently smooth parameter- dependent center manifold M cu (e) of the normally hyperbolic relative 
equilibrium SO(3)uq such that M{e) C M cu (e) and M cu (e) is diffeomorphic to 50(3) x V* for 
e > small. For e = 0, M cu (e) is SO (3) -invariant, and for e > small, M cu (e) is SO (2)- 
invariant. Furthermore, there exist sufficiently smooth functions f\ : V* — ► so(3), fa: V* — > V*, 
Gpert '■ SO (3) x Vif x [0, oo) — > Vif and a sufficiently smooth vector field G x peTt : SO (3) x V* x [0, oo) — ► 
T(50(3)) with Gp ert (.,q,e) being SO(2)-equivariant, G^ ert (.,q,e) being SO (2) -invariant such that 
any solution of 

A = Af x {q) + eG l pert {A,q,e) 

q = f2(q)+sG 2 pert (A,q,e), { ' ' 

on 50(3) x V* corresponds to a solution of the reaction- diffusion system \2.1\i on M cu (e) under 
the diffeomorphic identification for e > small. Also, /i(0) = Xq, /2(0) = and a(D q f 2(0)) = 
cr{Q*L\v t ), where Q* is the projection onto V* along T Uq (SO{3)uq). Moreover, is a hyperbolic 
equilibrium in V* for the second differential equation in \3. 1]) for e = 0. 

Henceforth, we will identify M(e) with SO{3) and M cu (e) with 50(3) x V* for e > small. 
Therefore, we will talk about the semifiow on 50(3) and on 50(3) x V* for e > small. 

Proposition 3.3. XTie general form of the SO(2)-equivariant perturbation on M cu (e) given by 
Theorem V3.2\ is 

G l p ert {A,q,e) = Ak(A,q,e), 
Gl ert (A,q,e) = h(A,q,e), 

where A G 50(3), q G V*, e > is small, and k: 50(3) X F* X [0,e ) -> so(3) and /i: 50(3) x7,X 
[0,£o) — ► so(3) are sufficiently smooth SO (2) -invariant functions for e G [0, £0) anc ^ 9 G V* raft 
eo > small. 

The proof of Proposition 13.31 is presented in Section H3 Generically, M(e) corresponds to a 
hyperbolic equilibrium q e G F near for e > small. Therefore, we can restrict the study to the 
50(2)-equivariant dynamics on 50(3), by substituting q = q £ in the differential equations (|3.1j) . 

Theorem 3.4. 77ie perturbed differential equations on M(e) given by Theorem V3. 2\ are 

A = A[X + eg(A,e)] (3.3) 

where A G 50(3), e > is small and, up to a constant matrix in so(3) that depends sufficiently 
smoothly on e, g(A, e) = k(A, q £ ,e) for k defined in Provosition \3.3l Thus g(., e) is SO (2) -invariant 
for e > small. 

The proof of Theorem 13.41 is presented in Sectional 
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4 The Analysis of the Projected Dynamics on the Orbit Space 

SO(3)/SO(2) 

Let x\ and x\ be the intersections of the line with direction the vector Xq with the unit sphere S 2 ; 
x i = anc ^ x 2 = ~]X^\Xq. Let and x® be the intersections of the line with direction the 

vector Q with the unit sphere S 2 ; = Q and x® = — Q . 
We consider the left action 9: SO (3) x SO (3) -> 50(3) given by 

0(^,O) = O^ 1 . (4.1) 

We substitute C = A' 1 in (ESJ). Since 3 = j 3 = AO = AC + AC, it follows that C = -CAC. 
Therefore, the differential equations (|3.3j) are equivalent to the following differential equations 

C = -[X Q + eg(C,e)]C, (4.2) 

where O G 50(3), e > small and 5 is a sufficiently smooth function such that, for e > small, 
g(C,e) = g(C ,e) is invariant under the restriction of the action 9 to the subgroup 50(2) of 
50(3), with g defined in Theorem 13.41 

Let us denote the flow associated to the differential equations (|4.2[) by It is clear that *f?(t,C,s) = 
C ,e)] . The differential equations (|4.2[) have the following property: 

• for s = 0, they are 50 (3)-equivariant under the action 6; 

• for e > small, they are 50 (2)-equivariant under the restriction of the action 9 to the 
subgroup 50 (2) of 50(3). 

We recall the notion of orbit space of 50(3) under the action 9 (see 0) and some of its properties 
that we will use to project the differential equations (|4.2|) from the phase space 50(3) onto the 
orbit space of 50(3) under the action 9. 

Definition 4.1. The quotient of SO (3) under the equivalence relation A ~ B iff A and B are on 
the same group orbit under the action 9 defined in \4-l\) and restricted to 50(2) is called the orbit 
space of SO (3) under the action 9 and it is denoted by 50(3)/50(2). It is the left coset of SO (3) 
modulo SO (2) under the action \4-l\j - 

Theorem 4.2 ( 123 EHJ ) . There is a unique C°° manifold structure for SO (3)/ S 0(2) such that 
the projection tt: SO (3) — > 50(3)/50(2) defined by vr(O) = C ■ SO (2) is C°° and such that there 
exist local smooth sections of 50(3)/50(2) in 50(3); that is, for any A ■ SO (2) £ S0(3)/ SO(2) 
there exists a local smooth function tt s from a neighborhood U of A ■ 50(2) into 50(3) such that 
tt o tt s = id on U. 

Theorem 4.3 ( j3 EHl |29j). The orbit space of SO (3) under the action 9 is diffeomorphic as a 
manifold to the unit sphere S 2 . The diffeomorphism can be chosen such that 

(5: SO(3)/SO(2) -» S 2 ,f3(A ■ 50(2)) = Ax® for any A G 50(3). (4.3) 

Definition 4.4. The projected flow on 50(3)/50(2) associated to the flow is defined by 

*(t, ?r(C), e) = 7r(*(t, C, e)) for any e > small , C G 50(3) and t > 0. (4.4) 
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The projected flow on 50(3)/50(2), 1 J r , is well-defined and sufficiently smooth. From [2U] and 
j2Zj) we know that every sufficiently smooth vector field on 50(3)/50(2) lifts to a sufficiently 
smooth 50 (2)-equivariant vector field on 50(3), unique up to vector fields tangent to the 50(2)- 
orbits. 

Definition 4.5. Using the diffeomorphism (3 defined in j4-3\) , we define the following flow on S 2 ; 

Cxf, e) = /9(*(t, ir(C), e)) for anye>0 small , C G 50(3) and t > 0. (4.5) 
The flows \I/ and are topologically equivalent. 

Theorem 4.6. The differential equations given by j4-2\ ) projected on the orbit space 50(3)/50(2) 
have the following form (after taking into account Theorem \4-3\ and Definition ^. 5| ): 

x = -[X + sg s (x,e)]x, (4.6) 

where x G S 2 and e >0 is small; that is, the vector field given by \4-b\ ) is the vector field associated 
to the flow ^i(t, Cx® , e). The function g s (x, e) is given by 

g s (x,e) = g(C,e), where x = Cx® G S 2 ,e > is small andg is defined in \4-2\ . (4.7) 
The proof of Theorem 14.61 is presented in Section H3 

Remark 4.7. 1. For e = 0, the differential equations J^.6| ) have reflectional symmetry with 
respect to the equatorial plane orthogonal to (the line containing) the vector Xq. 

2. For e > 0, generically, the differential equations \4-b\ ) have no reflectional symmetry with 
respect to the equatorial plane orthogonal to (the line containing) the vector Jo- 
Proposition 4.8 (e = 0). 1. For e = the solutions of the differential equations \3. 3\) are of 
the form 

A(t) = A e Xot , where Aq G 50(3). (4.8) 

2. For e = the solutions of the differential equations \4-b\ ) are the equilibrium points x® and x®, 
and the jj^-periodic solutions of the form x(t,xo) = e~ Xot XQ, whenever xq / x\ and xq ^ x\. 

3. The correspondence between the solutions of the differential equations h3.3\) and the differential 
equations j4-b\ ) is as follows: 

(a) The solution A(t) of \3. 3]) corresponding to Aq such that xf = Aqx\ projects onto the 
equilibrium x\. 

(b) The solution A(t) of 13.3)) corresponding to Aq such that xf = Aqx® projects onto the 
equilibrium x®. 

(c) The other solutions of i3.3\) with the initial condition Aq projects onto the periodic so- 
lutions of l4-b\ ) having the initial condition xq such that x® = AqXq. 

The proof of Proposition 14.81 is presented in Section HO For e = 0, the flow given by the 
differential equations ()4.2|) is the only 50 (3)-equivariant flow that can be obtained when we lift 
from the differential equations gSJ) on 50(3)/50(2) to 50(3). 
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line containing Q 



Figure 4.1: The Sphere 52 



Proof. The differential equations (|4.6|) for e = are x = —Xqx and its flow is given by ^(t,C ■ 
5O(2),0) = e~ Xot C ■ SO{2). We get vr(*(t,C,0)) = ¥(i,C • 50(2), 0) = e-*°*Cj 50(2). It 
follows that \l/(i,0, 0) = e _ ^- *Ce T ^' C '^. This is 50 (3)-equivariant under the action 9 if and only 
if r(t,C) = for t > and C G 50(3); that is, *(i,C,0) is the 50 (3)-equivariant flow given by 
differential equations (|4.2j) for e = 0. □ 

A geometrical interpretation of the dynamics on M(0) = 50(3)iio induced by the reaction- 
diffusion system ()2.1j) is given below. 

Proposition 4.9 (e = 0). We will consider two distinct unit spheres S\ and 52- Let the points 
and x® be on sphere S2 as they were defined at the beginning of this section. Also, let the points 
x\ and xS] be on sphere S\ as they were defined at the beginning of this section. For e = 0, the 
semiflow $ on SO(3)uo contains only rotating waves. We have: 

1. The rotating waves for the reaction- diffusion system \2. 1\) on SO{2>)uq with the unit frequency 
vectors X on a circle around Q of the unit sphere S2 project onto the periodic solution of 
i4-6\ )> that has as an orbit the circle around Xq of the unit sphere S\ passing through the 
point xq £ Si such that cos(Z(X ,Q)) = cos(Z(Xq,xo)). 

2. The rotating waves for the reaction- diffusion system h2.1\) on SO(3)uq with the frequency 
vector \Xq\ Q project onto x®. The rotating waves for the reaction- diffusion system \2.1\) on 
50(3)uo with the frequency vector — \X$\ Q project onto x^- 

3. Let C\ and C2 be two distinct circles on the unit sphere S2 with the centers on the line having 
the direction vector Q and such that d(0, C\) = d(0, 02). Let us denote by C[ the circle on 
S\ onto which project the rotating waves on SO(3)uo with the unit frequency vectors on C\. 
Let us denote by C 2 the circle on 5i onto which project the rotating waves on SO(3)uq with 
the unit frequency vectors on Ci- Then d(0, = d(0,C 2 ) and C x 7^ C 2 , where d is the 
Euclidean distance in R 3 . 

The following two figures illustrate Theorem 14.91 The Sphere 52 (Figure (|4.1|) 1 and The Sphere 
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the line containing XO 



Figure 4.2: The Sphere S\ 



Si (Figure (|4.2|l ). The proof of Theorem 14.91 is presented in Sectional 

Let B G 50(3) such that Bx\ = (0,0, 1) T . Let g ss {s,e) = Bg s (B- 1 s,e)B- 1 , where g s is defined 
in (|4.7|) and for any s = [x\, x 2 , x 3 ) G S 2 and e > small, 

/ -Gi(xi,x 2 ,x 3 ,e) G 2 (x 1 ,x 2 ,x 3 ,e) 

g ss (s,e) = g ss (xi,x 2 ,x 3 ,e) = (n(xi, x 2 , x 3 , e) -G 3 (xi, x 2 , x 3 , e) 

\ -G 2 (xi,x 2 ,x 3 ,e) G 3 (xi,x 2 ,x 3 ,e) 

Let 

a(4>, 9, e) = Gi(sin cf> cos 6, sin^sin^, cos cf>, e), 
b((j),6,e) = G 2 (sm(j)cos9,sin(psm6,cos(j),e), 

and 

c(4>,8,e) = G 3 (sin 4> cos sin<f> sin cos e), 

where G [0,27r), </> G [0,7r], e > small. The analysis of the perturbed differential equations (|4.6|) 
on the orbit space SO(3)/ SO(2) for e > small gives us the first main result of this article: 

Theorem 4.10 (e > 0). 1. There exist two sufficiently smooth branches of equilibria for \4-b\ ) 
defined for e > small, one branch x (e), near x® and the other one x 2 (e), near x\, such 
that x l (0) = and x 2 (0) = x%. 

Generically, we have ^^-(0,0,1,0) — ^^-(0,0,1,0) 7^ 0. The stability of the generically hy- 
perbolic equilibria x 1 (e) is as follows: 

(a) if ^^-(0,0,1,0) — ^^(0,0,1,0) < 0, then x 1 (e) is locally asymptotically stable for the 
differential equations j4-6\ ); 

(b) if 7%^{0, 0, 1, 0) — 7^(0, 0, 1, 0) > 0, then x 1 (e) is unstable for the differential equations 

Generically, we have ^^-(0,0,-1,0) — ^=^(0, 0, — 1, 0) 7^ 0. The stability of the generically 
hyperbolic equilibria x 2 (e) is as follows: 
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(a) if ^-(0,0, —1,0) - ^-(0, 0,-1,0) < 0, then x 2 (e) is locally asymptotically stable for 
the differential equations \4-b\) ; 

(b) if fgKO, 0, -1, 0) - fjKO, 0, -1, 0) > 0, then x 2 {e) is unstable for the differential equa- 
tions \4-b\ l- 



2. Let 

f27T 



1(4,) = / [b((f),t,0) cost - c{4>,t, 0) sin t] dt 
Jo 

for (p £ (0,ir), and consider 

(sin cos 8q \ 
sin 0o sin I with 0o £ (0, tv), 9q S [0, 2-k). For small e > 0, i/ie periodic solu- 
cos 00 / 

£ion x(t,xo) = e~ Xot xo of the differential equations \4-b\ l obtained for e = is continuously 
deformed when e increases to another periodic solution x £ (t) with period T(e) = + 0(e) 
if 

J(0 O ) = and j'(0 o ) + 0. 

AZso x°(i) = x(t, Xq). 

The stability of the periodic solution x £ (t) is as follows: 

(a) if t (4>o) < 0, then x £ (t) is locally asymptotically stable for the differential equations 

(b) if i' (0o) > 0, then x £ (t) is unstable for the differential equations \4-bty . 

Using the Poincare-Bendixson theorem on a sphere ([18]), we have a complete picture of the phase 
portrait for the differential equations j4.6\ ). 

The matrix B £ SO(3) such that Bx\ = (0,0, 1) T is not unique, but it can be shown that 
Theorem 14.101 does not depend on the choice of the matrix B. Therefore, in the statement and 
proof of Theorem I4.1U1 without loss of generality we may assume that = (0,0, 1) T . The proof 
of Theorem I4.1UI is presented in Section EJ 

Let e > be fixed and <fto 6 (0, it) such that it corresponds to the -p^-periodic solution x(t, xq) = 

e~ ot XQ of the differential equations (|4.6|) for e = 0. Then, the condition I(0q) = is necessary for 
the persistence of p^j-periodic solution x(t,xo) to T(e)-periodic solution x £ (t). 

5 The Analysis of the 5'0(2)-Equivariant Dynamics on the Phase 
Space 50(3) 

We translate the results of Section 01 which are valid for the orbit space SO(3)/SO(2), to the phase 
space 50(3) (Proposition 15. lj) or Y a (Theorem 15. 2 j) . 

Proposition 5.1 (e > 0). For £ > small we have: 

1. The equilibrium x x (e) of the differential equations \4-o\) gives the SO(2)-orbit of a solution 
of the differential equations \3. having the form A(t,£) = e ai ( £ ^*A £) where xf = A^ 1 ^) 
and ai(e) = |Xo| +(9(e). The equilibrium x (e) is generically hyperbolic, therefore the corre- 
sponding SO(2)-orbit is generically normally hyperbolic; 
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2. The equilibrium x 2 (e) of the differential equations \4-b\) gives the SO(2)-orbit of a solution 
of the differential equations \3.3\) having the form A(t,e) = e a2 ( £ ^*A e) where xf = A £ x 2 (e) 
and 02(e) = — \Xq\ + 0(e). The equilibrium x 2 (e) is generically hyperbolic, therefore the 
corresponding SO(2)-orbit is generically normally hyperbolic; 

3. If I(0o ) = and I (<po) ^ 0, where I and 4>q are defined in Theorem \4-l(A then any T(e)- 
periodic solution x £ (t) of j4.b\ ) gives the SO(2)-orbit of a solution of the differential equations 
l/Oj) having the form A(t,e) = e^ Qt B(t,e), where /3(e) = 0(e) and B(t,e) is a T(e) = 

+ O(e)- periodic function such that xf = B(0,e)x e (0). 

The orbital stability of the above solutions is the same as the stability of the solutions of j4.6\ ) from 
which they were obtained (see [11). Using the Poincare-Bendixson theorem on a sphere and 
0/, we have a complete picture of the phase portrait for the differential equations \3. 3\) . 

The proof of Proposition 15.11 is presented in Section |H1 
The second main result of this article is the following theorem: 

Theorem 5.2 (The Effects of Forced Symmetry-Breaking). For small e > we have: 

1. The SO(2)-orbit of the rotating wave 3>(t, Aquq,0) = e^ x °^Aouo is deformed as e increases 
to the SO(2)-orbit of a rotating wave having as the frequency vector (\Xq\ + 0(e)) Q , where 
Q = AqXq; 

2. The SO(2)-orbit of the rotating wave $>(t, AqUo,0) = e~^ x °^AQUQ is deformed as e increases 
to a rotating wave having as the frequency vector (— \Xq\ + 0(e))Q ; 

3. The SO(2)-orbit of the rotating waves &(t, AqUq,0) = e xt A()Uo having the frequency vectors 

X = AqXqAq 1 different from \Xq\ Q and — \Xq\ Q , and satisfying I((fto) = and I (cfio) 7^ ; 
where I and <j)Q are defined in Theorem \4-l(J[ transform as e increases into SO(2)-orbits of 
modulated rotating waves having primary frequency vectors of the form 0(e) Q and the relative 
period of the form T(e) = pj^r + 0(e). 

The orbital stability of the above solutions of the reaction- diffusion system \2.1\) is related to the 
orbital stability of the rotating wave *&(t,uo,0) and to the stability of the solutions of \4-b\l from 
which they were obtained as follows: 

• If the rotating wave <&(t,uo,0) is orbitally unstable, then the above solutions of the reaction- 
diffusion system i2.1\) are orbitally unstable; 

• // the rotating wave uq, 0) is orbitally stable (that is E cu = T uq (SO(2>)uq) ), then the above 
solutions of the reaction- diffusion system \2.1\) have the same orbital stability as the stability 
of the solutions of \4-b\l from which they were obtained. 

The proof of Theorem 15 .21 is presented in Section H3 In 0, we present some examples that may 
explain similar phenomena to those observed for planar spiral waves (boundary drifting, anchoring 
at the localized inhomogeneities, repelling from the localized inhomogeneities) , and another new 
phenomenon that appears due to the curvature of the sphere (anchoring at the opposite site). 
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6 Proofs of Theorems 

Proof of Proposition^^ There is nothing to prove for Gp ert (A, q, e). 
Let k{A,q,e) = A~ l Gl, ert (A,q,e) for A G 50(3), q G V* and e > small. 

Since G l pert {A, q, e) G 7^(50(3)) = A - so(3), then fc(A,g,e) G so(3) for any A G 50(3), q G V* and 
e > small. Using the 50 (2)-equi variance of G pert (.,q,e), we get that the function k(.,q,e) is 
50 (2)-invariant for e > small and q G V*, that is 

fe(^i,? lE ) = (e^A^G^e^A^e) = A^e-^e^G^fi,^) = A^^g.e) = fc(^,g,e) 
for any A G 50(3), q€V+,0€ [0, 2vr) and e > small. □ 

Proof of Theorem \'J.4\ For e > small, M(e) corresponds to the hyperbolic equilibrium q £ G U* 
of the second differential equation of (|3.1|) . Therefore, we only consider the 50 (2)-equi variant 
dynamics on 50(3) obtained by substituting q = q e in the differential equations (|3.1|) . Using 
Proposition 13.31 we get 

A = A[/ 1 ( fe )+eA:(A &,e)]. (6.1) 

Since q e = 0(e) and f\ is sufficiently smooth, we get fi(q e ) = /i(0) + eH(e) = Xo + eH(e) with 
ff(e) G so(3). The differential equations (|6.1|) become 

A = A[X + eH{e)+ek(A,q £ ,e)}. (6.2) 

Let g(^4, e) = H(e) + fe(A, q £ ,e). We get the differential equations ()3.3|) . where g(-,e) is 50(2)- 
invariant, since k(.,q e ,e) is 50(2)-invariant. □ 



x 



Proof of Theorem \4-o\ It is enough to prove that there exists a sufficiently smooth function g s : S 
[0, £0) — > so(3) for £0 small such that 

«f(C e) = g S (Cx 1 Q , e) for any O G 50(3) and e > small . (6.3) 

Since S 2 is compact, it is enough to prove ()6.3|) locally on S 2 . We may then use the partition of 
unity to get a sufficiently smooth function g s globally defined on S 2 . 

Let Cq G 50(3) be fixed, but arbitrary. There exists a smooth local section ir s of 50(3)/50(2) in 
50 (3) near C ■ SO (2) (see Theorem E2J) . We have: 



7T O 7T S = id 

7r{7r s {C ■ 50(2))) = C ■ 50(2) 
=> 7T S {C ■ 50(2))50(2) = C ■ 50(2) 

C = ir s (C ■ 50(2)) • C*{C) for C*(C) G 50(2) near O • 50(2). 



(6.4) 



Let us define the function 

g s = g o 7r s o near Oo^q. (6-5) 
Then g s is a sufficiently smooth function defined locally near CqXq on S 2 . Also, 

g s (Cx 1 Q ,e)=g(w s (C-SO(2)),s)=g(7r s (C-SO(2))-C*(C),e)=g(C,s) (6.6) 

for C near Co and e > small, where for the second equality we use the invariance of g(-,e) under 

the restriction of the action 6 to 50(2) and for the third equality we use (|6.4j) . 

The function g s is well-defined since the function g is SO (2)-invariant under the action 9 restricted 
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to 50(2); that is, for any 9 G [0,2vr), C G 50(3) and e > small, we have g s {Ce Q0 x l Q , e) = 
g(CeQ e ,e)=g(C,e)=g s (Cx 1 Q ,e). 

Also, taking into account that the action of 50(3) on S 2 is transitive, we can look at the function 
g s as being defined on S 2 . 

We check that the vector field given by the differential equations IJ4.fi|) is associated to the semiflow 
*i. We have Wx(t,Cx l Q ,e) = (/? o 7r)(¥(i, C, e)) = ^(t,C,e)x 1 Q . Therefore, 

-^±(t,Cx l Q ,e) = -m,C,s)x 1 Q ) = - w (t,C,s)x l Q 

= -[X + eg(^(t,C,e),e)]^(t,C,e)x 1 Q (g7) 
= -[X + eg s (*(t, C, e)x 1 Q , e)] C, e)x l Q 
= - X + eg s (¥ 1 (t,Cx 1 Q ,e),£^ *i(t,Cx^e) 

for any C G 50(3), t G [0, oo) and e > small. □ 



Proof of Proposition ^. 8\ Conclusions (1) and (2) are obvious. The third conclusion results from 
the definition of the flows * and *i, as well as C, e) = [<&(t, C' 1 , e)]^ 1 . □ 

Proof of Theorem \4.y\ Conclusions (1) and (2) result from Theorem 13.21 and Proposition 14.81 

Using Theorem 13.21 we get that the rotating waves for the reaction-diffusion system (|2.1j) for e = 

are of the form &(t,A u ,0) = A e Xot u for A G 50(3) and t G [0,oo). 

> 

Therefore, the primary frequency vector of $>(t, Aoiio,0) is ^o^o^o 1 - By Proposition 14.81 the 
rotating wave <&(i, AqUq,0) projects onto xq G S 2 such that x® = Aqxq. 

Let X = j^A X A 1 . Then »$ = p^y^o, we get A x? = X. Also, (xf ) T = x%A%. Therefore, 
cos(Z(X,Q)) = (x®) T X = xlAlA x\ = a#s° = cos(Z(x , X )). 

The third conclusion results from the first conclusion. □ 

Proof of Theorem \4-l(J\ For e = 0, the differential equation (|4,6|) has two equilibria, x® and x\. 
Using the Implicit Function Theorem , we will prove the persistence of these two equilibria for 
small e > 0. 

We have defined B G 50(3) by Bx\ = (0,0, l) T . Then Bx\ = (0,0, -1) T . Let s = Bx. Then, 

s = Bx = - [BXqB- 1 + eBg s {B- l s, e)B-~ l ] s, (6.8) 
and, if we write Xoo = BXqB^ 1 and g ss (s,e) = Bg s (B~ l s,e)B~ 1 , then we get 

s = -[X 00 + eg ss (s,e)]s, (6.9) 

Since BXq = \X$\ (0,0, 1) T , we have 

-^oo = BXqB~ 1 = \Xq 

Let us define the function 

F(s,e) = - [Xqo + eg ss (s,e)] s, where F : S 2 x [0,e ) TS 2 for e > small . (6.10) 
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It is clear that F((0,0, 1),0) = 0. We will compute (A^^oai)^) : T (o,o,i) s2 ~* TS2 - r2 - 
Let s = (xi,x 2 ,x 3 ) G S 2 and recall that 



g ss (s,e) = g bb {x 1 ,x 2 ,x 3 ,e) 



ss, 



-G 1 (x 1 ,x 2 ,x 3 ,e) G 2 (xi,x 2 ,x 3 ,e) 

Gi(xi,x 2 ,x 3 ,e) -G 3 (xi,x 2 ,x 3 ,e) 

-G 2 (xi,x 2 ,x 3 ,e) G 3 (xi,x 2 ,x 3 ,e) 



For (xi,x 2 ,x 3 ) E S 2 near (0,0,1), one has x 3 = \J\ — x\ — x\ > for ||(xi,x 2 )|| > small enough. 
If we substitute this expression of x 3 into the differential equations ()6.9|) . we have: 



x\ \ 






-\X \ 





x 2 


1- 


.( i l 








x l x 2 / 













+e Gi(xi,x 2 , y/l -x\- s|,e) 



-G 2 (x 1 ,x 2 , y/l 



-G 1 (x 1 ,x 2 , y/\-x\- s|,e) G2(xi,X2,\/ l - xf - x|, e) 
-G 3 (xi,x 2 , v 7 ! - xf - x|,e) 





;,e) G 3 (xi,x 2 , V 1 ~ ^i ~ x i> e ) 



X2 



or 



.r-i 



|-^o I ^2 +eGi(x,e)x 2 - eG 2 (x,e)v / l - x? 
- |X | xi - eGi(x 1 e)xi + eG 3 (x,e)y / T 



-2> 



x 2 > j 



(6.11) 



(6.12) 



eG 2 (x,e)x\ - eG 3 (x, e)x 2 , 



where x = {x\,x 2 , y/l — x\ — x|) . 

If the first two equations in the system (|6.12|) are satisfied, then the third one is immediately 
satisfied. Therefore, we get the system 



x\ = Hi(x 1 ,x 2 ,e), 
x 2 = H 2 (xi,x 2 ,e), 



(6.13) 



where 



H 1 (x 1 ,x 2 ,e) = \X \ x 2 + e 



Gi(xi,x 2 , Jl-xf- x\,e)x 2 



G 2 (xi,x 2 , Jl — x\ — x 2 ,e)Jl — x\ — x 



(6.14) 



and 



H 2 (xi,x 2 ,e) 



\X \ xi -e 



Gi(xi,x 2 , Jl-xf- x%,e)xi 



G 3 (xi,x 2 ,Jl -x[- x\,e)J\ -xf-x 



(6.15) 



Let H(xi,x 2 , e) = (Hi(xi,x 2 ,e), H 2 (xi,x 2 ,e)). We have H(0, 0, 0) = (0, 0) and the linearization of 
H(xi,x 2 ,e) about (0,0) at e = is 



D (XUX2) H(0,0,0) 



\X \ 
\X \ 
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which is an invertible matrix. Using the implicit function theorem there exists a sufficiently smooth 
branch (x\(e), a4( e )) near (0,0) such that (x\(0),xl(0)) = (0,0), and Hi(x\(e),xl(e),e) = and 
H2(x\(e), x\{e), e) = for e > small. 
We have 

x\(e) = x\\E + 0(e 2 ) and x\{e) = x 2 2^ + 0(e 2 ), for small e > 0. 



By using implicit differentiation we get 



8H\ dx\ _|_ 



dH\ dx2 



dx± de 8x2 de 
8H2 dx\ _|_ 8H2 dx2 
dxi de 8x2 de 



dHi 
de ' 

8H 2 
de ■ 



Generically, we have 



G 2 (0, 0, 1, 0) / and G 3 (0, 0, 1, 0) / 0. 
From (|6~T7|) it follows that at e = we have 



(6.16) 

(6.17) 
(6.18) 



and 



dx 2 
de 



am 

de 
dH 2 
~5T 


am 

dx2 

dm 

dx2 


dm 


dm 


dxi 


dx2 


dm 


dm 


dxi 


dx2 


dm 


dm 


dxi 


de 


dm 


dm 


dxi 


de 


dm 

dxi 

dm 

dxi 


dm 
|E 

8x2 



(6.19) 



(6.20) 



Since §f(0,0,0) = §g (0,0,0) = 0, |g(0,0,0) = -fff (0,0,0) = \X \, ^(0,0,0) = -G 2 (0, 0,1,0) 
and 2^(0, 0, 0) = G 3 (0, 0, 1, 0), it follows that 



xu 



^G 3 (0, 0,1,0), 



\x 



X22 = ^02(0,0,1,0). 



So, for e > small, we get 
x\{e) 



-r^G 3 (0, 0, 1, 0)e + 0(e 2 ), x\{e) = 
A o l A o| 



G 2 (0,0,l,0)e + O(e 2 



(6.21) 



(6.22) 



Therefore, there exists a sufficiently smooth branch s (e) = (x\(e) , x^e) , xl(e)) near (0,0,1) of 
equilibria of the differential equations Q6.9JI . such that 

s i(0) = (x\(0),xl(0),xl(0)) = (0,0,1) and F{s 1 (e),e) = (0,0,0) for e > small. 
For e > small, there is a sufficiently smooth branch x x (e) G S 2 of equilibria for the differential 
equations (|4.6|) such that x x (0) 



'V 



In the same way, for e > small there exists a sufficiently smooth branch x 2 (e) of equilibria for the 
differential equations l)4.6j) such that x 2 (0) = x®- 

The stability of the equilibrium x 1 (e) is the same as the stability of the equilibrium {x\(e) , x\(e)) for 
e > small, and the stability of the equilibrium x 2 (s) is the same as the stability of the equilibrium 
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(x 2 (e), x|(e)) for e > small. 

Let us now obtain the stability of the equilibrium (x\(s), x\{e)). 

For e > small, the linearization of H(x%, X2, e) around (xj(e), x^e)) is given by 

D,„ m) H(,\(s)y^U) = ( _ 0(£) ™ ;°< £ > ) , (6.23) 



where 



and 



dH 

a u ( £ ) = -^(x\{e), x\{e), e) = a u e + 0{e 2 ) since a u (0) = (6.24) 
dH 

= -r. — (x\(e) , x\{e) , e) = a 22 e + 0(e 2 ) since a 22 (0) = 0. (6.25) 



<9x 2 

The eigenvalues of the linearization of .ff(xi, x 2 , e) around (x\(e), x\{e)) satisfies 



on(e)-A \X \+O(e) 
-|*o|+0(e) a 22 (e)-A 



0. (6.26) 



Then, it follows that (a u (e) - A) (022(e) - A) - (\X \ + 0(e))(- \X Q \ + O(e) = or 
A 2 — (on(e) + a 22 (e))A + \X \ 2 + 0(e) = 0. We get the eigenvalues 



aii(e) + a 22 (e) ... 1 , 
Ai, 2 (e) = o ± Hl^ol + 0(e)) 



an + a22 „, 2 ^ 
o e + 0(e ) 



±i(|X |+O(e)). (6.27) 



Generically, an + a 22 7^ 0. In fact, (x\{e),x\{e)) is generically a hyperbolic equilibrium for the 
differential equations (|4.6|) . If an + a 22 < 0, then (x\(e), x\{e)) is an asymptotically stable equi- 
librium for the differential equations (|4,6|) . and if an + a 22 > 0, then (x\(s), x\{e)) is an unstable 
equilibrium for the differential equations (|4.6j) . 

We express an and a 22 in terms of Gi((0, 0, 1), 0), G 2 ((0, 0, 1), 0) and Gs((Q, 0, 1), 0) and their par- 
tial derivatives at (0,0,1). 

The Taylor expansion of Hi(x\, X2, e) about (0,0,0) up to terms of order 2 (h.o.t. contains all the 
terms of order > 3 in xi, X2, e) is: 



Hi{xi,X2,e) = al + ajxi + a 2 x 2 + a\x\ + a^x 2 , + a\xix 2 
+ ag£ + a\x\£ + a\x2£ + ag£ 2 + h.o.t. 



(6.28) 



Since aj = 0, a\ = 0, 03 = 0, a\ = and a\ = 0, we get 

H\(xi,X2,£) = a\x2 + a\e + a\x\e + a§x 2 e + a^e 2 + n.o.t. (6.29) 
The Taylor expansion of £f 2 (xi, x 2 , e) about (0,0,0) up to terms of order 2 is: 

H2(x\,X2,s) = a 2 , + a 2 xi + a 2 x 2 + a 2 x 2 + a|x 2 + a 2 xix 2 

22 2 2 2 7 ^6.o0j 
+ a 6 e + a 7 xie + a 8 x 2 e + a 9 e + h.o.t. 

Since a 2 , = 0, a 2 = 0, a 2 = 0, a 2 = and a 2 = 0, we get 

i/ 2 (xi, x 2 , e) = a 2 xi + ag£ + a 2 x\£ + agX 2 e + a^e 2 + h.o.t. (6.31) 
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We have 



and 



We have that 



an 



by the definition of Hi and 



022 



dHi 

dxi 

dH 2 

dx 2 



(x\{e),xl(e),e) = a\e + o(e 2 ) 



(x\(£),xl(e),e) = a 2 8 e + o(e 2 ). 



d 2 Hi 

dxide 



(0,0,0) 



dG 2 
dxi 



(0,0,1,0), 



d 2 H 2 



(0,0,0) 



0G, 



(0,0,1,0), 



(6.32) 
(6.33) 
(6.34) 

(6.35) 



dx 2 de dx 2 
by the definition of H 2 . 

Generically, we have a\ + a| ^ 0. If a\ + a| < 0, then (xj(e), x\{e)) is an asymptotically stable 
equilibrium for the differential equations ()4.6|) . and if ay + a| > 0, then (x\(s) , x 2 (e)) is an unstable 
equilibrium for the differential equations ()4.6|) . 

In the same way, we can establish the stability of the equilibrium x 2 (e). This ends the proof of the 
first conclusion of Theorem 14. 101 

For e = the differential equation (|4.6|) has po _ periodic solutions of the form x(t, xq) = e~ Xot xo, 
where xq / x\ and xq / x 2 . 

Let us fix a point xo on the unit sphere such that xq / x® and xo 7^ x 2 . 

We will prove the persistence of the periodic solution x(t,xo) if some conditions are satisfied. 
Making the same change of variable as before, we get the system ()6.9|) and the periodic solution 
x(t, xo) becomes s(t, sq) = e~ x ° ot so, where so = Bxq and s$ / (0, 0, 1), sq / (0, 0, —1). 
For s 7^ (0, 0, 1) and s 7^ (0, 0, —1) we use the spherical polar coordinates 
sin <t> cos ( 

where (ft G (0,tt) an d 9 G [0, 2n). Recall that sq 



Let s 
and 



sm <p sm ( 
cos 6 



,9) on S 2 . 

sin cf)Q cos #0 
sin rf>o sin #0 
cos (bo 



for 6> G [0, 2vr) 



a(4>,9,e) = G 1 (sin cj) cos 9, sin cfi sin 9, cos <f>, e), 

b(4>,6,e) = G 2 (sin cf> cos 9, sin (p sin 9, cos (p,e), 

c((f>, 9,e) = G3 (sin cos 6*, sin </> sin cos <fi, e) 

G [0, 7t] 3 e > small and such that 

a(0,8,e) = Gi(0,0, 1, e), a(vr, 6», e) = Gi(0, 0, -1, e)for any 6> G [0,2vr), 

6(0, 0,e) = G 2 (0,0, l,e), 6(tt, 0,e) = G 2 (0, 0, -1, e)for any 9 G [0,2vr), 

c(0,9,e) = G 3 (0, 0,1, e),c(TT,9,e) = G 13 (0, 0,-1, e)for any 

The system ()6.9|) in (</>, 0) coordinates yields: 

^ (sin 4> cos 6>) 
^ (sin 4> sin 6>) 
|(cos0) 



+e 



G [0,2tt). 



- 


\*o 1 


l*o| 














—a 


{$,6,e) 





-6(0, e,e) 








^ sin 4> cos \ 


)} 


1 


sin <j) sin 






^ COS (j) / 
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> cos (ft cos 8 — 8 sin (ft sin 6 

> cos (ft sin 6 + 6 sin (ft cos 

— c/>cos 6 




- |X | -£a((ft,9,e) eb{(ft,6,e) 

\Xq\ + ea(4>,6,e) —ec((ft,d,e) 

-eb{4>,6,e) ec{(ft,8,e) 



</> cos <^> cos 6 — 6 sin (ft sin = [|Xo| + ea((ft, 6, e)] sin (ftsin — e6((/>, 0, e) cos 0, 
<^> cos (ft sin + # sin (ft cos 6* = — [|Xo| + ea((j>, 6, e)] sin^cos 6 + ec((ft, 6, e) cos (ft, (6.36) 
—(ft sin (ft = eb{(ft, 6, e) sin 0cos 6 — ec((ft, 6, e) sin (ft sind. 

We divide the third equation in (|6.36[) by sin (ft ^ 0, where (ft £ (0,vr) to get 

ift = -sb(cft, 6, e) cos 6 + ec(cft, 6, e) sin 6. (6.37) 

We substitute (ft given by (|6.37j) into the first two equations of (|6.36|) to get 

6 sin</>sin 6 = — \\X$\ + ea((ft, 6, e)] sind sin(ft + ec((ft, 6, e) cos 6 sind cos 

e) cos sin 2 6, 

6 sin (ftcos6 = — [\Xo\ + sa((ft,6,E)] cos 6 sin (ft + ec(0, 6, e) cos (ftcos 2 6 

+eb(tft,9,e)sin6cos6cos(ft. (6.38) 

If sin# ^ 0, we divide the first equation in (|6.38j) by sin</>sin#, to get 

6 = — [\Xo\ + ea((f), 9, e)] + eb(<ft>, 6, e) sin 6 cot (ft + ec((ft, 8,e) cos 6 cot (ft. (6.39) 

Thus, the second differential equation of (|6.38|) is satisfied. 

If sin 6 = 0, then cos 8^0 and we divide the second equation in (|6.38|) by sin (ft cos 6 to get 

= - [|X | + ea(</>, (9, e)] + eb((ft, 6, e) sin 8 cot (ft + ec((ft, 6, e) cos 6 cot (ft. (6.40) 

Thus, the first of the differential equations (|6.38j) is satisfied. 
Therefore, we obtain the system 



<ft = — eb((j), 6, e) cos 6 + ec{(ft, 8, e) sin 6, 

8 = — [\Xo\ + ea((ft, 8, e)] + eb(4>, 6, e) sin 6 cot <ft> + ec((ft, 8, e) cos 8 cot < 

or 

(ft = e[— b((ft, 6, e) cos 8 + c(4>, 8, e) sind] , 

6 = -\X \ + e{a((ft, 8, e) + [b(4>, 6, e) sind + c((ft, 8, e) cos 8} cot (ft}. 
For (ft 6 (0, 7r), 8 6 [0, 2ir) and for e > small enough, we can have that 

- \X \ + e{a((ft, 8, e) + [b{(ft, 6, e) sin 8 + c(<ft, 8, e) cos 0] cot </>} > 



3.41) 
3.42) 
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and we can rescale time, choosing r = r(t) such that 



We get 

Since 
and 

it results that 



dr 

— = - \x \ + e {a(0, 9, e) + [6(0, 9, e) sin 9 + c(0, 9, e) cos 9] cot 0} . (6.43) 

d<£ _ d4 dt_ _ e[—b(4>,B,e) cos 6+c(<t>,6,e) sin 6] 

dr ~~ dt dr ~ -\X \+s{a(4>,8,e)+[b(4>,8,e) sin 6+cl<j>,6,e) cos 6*1 cot 0} ' f« 44 N 

de _ dt _ 1 v u - I± ' i ; 

dr dt dr 

6(0,6>, £ ) =b((f>,9,0)+eh 1 (<l),e,e) 
c((p, 9, e) = c((p, 9, 0) + e/i 2 (0, 9, e), 



(6.45) 



f = p^r [6(0, 0, 0) cos 0-c(0,0,O) sin 0] + e 2 h(0, 0,e), 
de _ i 

dr ~~ L - 

Let us construct the Poincare map ( see |30| . that is the time 2n map) for the flow given by the 
differential equations 1)6.451) . The Poincare section associated to the flow given by the differential 
equations (|6.45|) is given by 9 = 0, G (0, n), with chosen independent of 9q G (0, tt). The 
Poincare map is given by 

P(0i,e) = 0(2tt,0 1 ,£), 



where 0(t, 0i, e) is the solution of the following initial problem after we relabel r = t, 

# = £ [6(0, t,0) cost - c(0,i,O) sin i] + e 2 /i(0, t, e), 
0(0) = 0i° 

Therefore, we get by using ()6.46|) 



(6.46) 



P(0i,e) =0(0,01, e)+ / -J(t,0i,e)dt 



2tt 



2tt 



'j H / <j t^t [6(0(t, 01, e), t, 0) cos t - c(0(t , 0i , e) , t, 0) sin t] 



+e 2 /i(0(t,0i,e),t,e)} (it 



or 



2 - E 



P(0i,e)=0i+ / — — [6(i, 0i,O) cost - c(t, 0i,O) sin t] dt + e h* (0i, e). 
7o l^o | 



We have that 



where we define 



P(0i,e) = 0i + — - /(0i) +£ 2 /i*(0i,£), 
l A o| 

/•27T 

1(0!) = / [6(0i, t,0) cos t-c(0i,t,O) sin t] dt. 
Then, 

P(0o,O)-0 o = O, 
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since x(t,xo) = e~ Xot XQ is a periodic solution of the differential equations (|4.6|) for e = 0, that 
is s(t,sn) = e~ x ° ot so is a periodic solution of the differential equations (|fi.9|) for e = and so = 
sin 0o cos $o 

sin 0n sin 8q ) . Let us consider the equation 



COS ( 



that is 
If 



P(0 1)£ ) - = TFT [/(0i) = 0. 



l(0i) +eh*(<f>i,e) = 0. 



I(0o) = and I (0 O ) ^ 0, 

then, using the implicit function theorem, we find for e > small a sufficiently smooth branch 0(e) 
of fixed points of P such that 0(0) = 0o- 
We have 

6(01 , t, 0) = G2 (sin 0i cos t, sin 0i sin i, cos 0i , 0) 

and 

c(0i , t , 0) = G3 (sin 0i cos i, sin 0i sin t, cos 0i , 0) . 



Therefore, 



db dGc2 

— — (0i, t, 0) = — — (sin 0i cos t, sin 0i sint, cos 0i, 0) cos 0i cos t 
dG 

+ — — (sin 0i cos t, sin 0i sin t, cos 0i , 0) cos 0i sin t 
0x2 

dG 

— - — (sin 0i cos t, sin 0i sin t, cos 0i , 0) sin 0i 

9x 3 



and 



dc dG 

— — (0i, t, 0) = — — (sin 0i cos t, sin 0i sint, cos 0i, 0) cos 0i cos t 

£701 £7Xl 

+ — — (sin 0i cos t, sin 0i sin t, cos 0i , 0) cos 0i sin t 

0x2 

dG 

— - — (sin 0i cos t, sin 0i sin t, cos 0i , 0) sin 0i . 
ox 3 

For e > small, we get a sufficiently smooth branch s £ (t) of periodic solutions for the differential 
equations (|6.9j) such that s°(t) = e~ x ° ot so. 

Thus, for e > small, we get a sufficiently smooth branch x e (t) of periodic solutions for the 
differential equations (|4.6|) such that x°(t) = e~ X()t xo. 

If / (0o) < 0, the fixed point 0(e) is locally asymptotically stable for P, and if / (0o) > 0, the fixed 
point 0(e) is unstable for P. 

The stability of the periodic solutions x £ (t) which persist for e > is the same as the stability of 
the fixed points 0(e) of the Poincare map P from which the periodic solutions x £ (t) are obtained. 
This ends the proof of the second conclusion of Theorem 14.101 □ 
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Proof of Provosition \5.1\ We prove the first conclusion for the equilibria x 1 (e). A similar proof can 
be done for the equilibria x 2 {e). 

Let A £ e SO (3) be such that xf = A e x 1 (e). We can choose a sufficiently smooth branch A £ , by 
taking C £ = (ir s o /J- 1 )^ 1 ^)) and then A £ = C £ x . We have C e xf = x 1 ^). 

Taking into account that ^(t,C,e) = [&(t, A, e)] _1 , with A = C _1 , we check that the differential 
equations (|4.2[) have a solution of the form 

C RW {t,e) = C £ e' ai(£)Qt with ai (e) = \X \ +0(e). 
By the definitions of the projected flows and we have C £ x® ,e) (e), then 

*(*, C £ ■ SO(2), e) = 7r(*(i, C e , e)) = C e • 50(2), 

which implies 



C7 w (t,e) = ¥(t,C e ,e) = C £ e T ^ e)Q , 



(6.47) 



where r(0, e) = 0. Without loss of generality, we assume that r(t,e) sufficiently smooth with 
respect to t. It follows that C (t, e) verifies the differential equations ()4.2jl . 
If we substitute Q6.47JI into the differential equations (|4.2j) . we get 



C £ ^(t,e)Qe^Q = - [x + eg(C £ e^^,e) 



Q e r{t,e)Q 



or using the S'0(2)-invariance of g(-,e), we get 



dr 



-[X Q + eg(C e ,e)]Q 

8t 
Ot 
dr 



(t,e)Q = -C- 1 [X +eg(C £ ,e)]C £ 



dt 



(t,e)Q 



C- L [X + eg(C £ ,e)]C £ 



or, using \Q\ = 1, we get 



dr . . 



\X + eg(C £ ,e)\ 
— (t,e)= constant = f Qi(e). 



(6.48) 



Therefore, 



r(t, e) = —a\(e)t + r(0, e) => r(t, e) = —a\{e)t. 
Also, we get a\{e)Q = C~ l [Xq + £g(C £ ,e)C £ ] and a%(s) is sufficiently smooth. Then, 

ai(e) = \X \ + O(e), 

since a\(0)Q = Cq 1 XqCq implies a\(0)Q = \Xq\ Q, namely ai(0) = \Xq\ because xf = C c j~ 1 x 1 (0), 
a?i(0) = x\. 

The stability issue for C RW (t, e) is proved in 0. It is now easy to get the results for the differential 
equations using the fact that *(t,C,e) = [*(£, ,4, e)] -1 , where ^ = C _1 . 
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This proves the first conclusion. Similarly, we can prove second conclusion of Proposition 15.11 
We now prove the conclusion (3) of Proposition 15.11 

Let D(t,e) G 50(3) be such that D(t,e)xf = x £ (t). Let T(e) = t^t + 0(e) be the period 
of the function x £ (t). We can choose a sufficiently smooth branch D(t,e), by taking D(t,e) = 
(ir s o (3^ 1 )(x £ (t)). Since x £ (t) is T(e)-periodic, we have D(t + T(e),e)x® = D(t,e)x®. Therefore, 
since the section ir s is local, we can only say that there exists a sufficiently smooth function g(t, e) 
such that 

D(t + T{s), e) = D(t, e)e^' £)Q . (6.49) 
Without loss of generality, we assume that g(0, e) = 0, otherwise we take D\(t, e) = D(t, e)e~ 9 ^°' £ ^ . 

Taking into the account that C, e) = [3>(t, A, e)]" 1 , where A = C _1 , we check that the differ- 
ential equations (|4,2[) have a solution of the form 

C MRW (t,e) = B*(t,e)e-^ t 

with /3(e) = 0(e), B*(0,e) = D(0,e) and B*(t,e) is T(e)-periodic. 

Let C\ = D(0,e). We construct B*(t,e) and (3(e). 

By the definitions of the projected flow and we have ^t(t,C^xf,£) = x £ (t) = D(t, e)x® or 
¥(i,C£ ■ SO(2),e) = ir(V(t,C^,e)) = D(t,e) ■ 50(2), which implies 

*(t, Cl,e) = D(t, e)e~ T(t ' £)Q with r(0, e) = 0. (6.50) 

Without loss of generality, we assume that r(t,e) is a sufficiently smooth function with respect to 
t. If we substitute (|6.50f) into the differential equations ()4.2(l . we get 

D = - [X Q + eg(D, e)} D + fDQ => fDQ = D~ x t> + D' 1 [X + eg(D, e)] D. 

Therefore, the function r(t, e) is sufficiently smooth in e and in t. We define 

C MRW (t,e) d ^ D(t,e) e -< t ^e^ t e-^\ where [3(e) = ■ 

Let us define B*(t,e) = D(t, e)e~ T ^ Q e^ Qt . Then, C MRW (t,e) = B*(t,e)e~^ Qt . 

We have that B*(T(e),e) = D(T(e), £ ) e -r(T(e),e)Q e p(e)T(e)Q _ D ( Q ^ £ yg(o,e)Q = D(0,e) = B*(0,e) 
by the definition of (3(e), relation (|6.49|) and the fact that g(0, e) = 0. 

Since r is sufficiently smooth in t and e and T(e) is sufficiently smooth, it follows from the definition 
of (3(e) that [3(e) is a sufficiently smooth function and [3(e) = + 0(e). 

We check that B*(t,e) is T(e)-periodic. If we substitute C = Be~ l3 ^ < ^ t into the differential 
equations (|4.2|) . we get 

Be -mQt + B(-P(e)Q)e-^ £)Qt = -[X + eg(B, e)] Be~^ Qt 



or 



B = (3(e)BQ- [X + eg(B, e)]B. (6.51) 



23 



Since C MRW {t, e) = B*(t, £ )e^^ Qt and C MRW (t, e) is a solution of the differential equations (|4~2"|) . 
we get that B*(t,e) is a solution of the differential equations (|6.51j) . Since B*(T(e),e) = B*(Q,e), 
we get that C*(t, e) = B*(t + T(e),e) is also a solution of the differential equations (|6,51|) such that 
C*(0,e) = B*(0,e). Therefore, the function B*(t,e) is T(e)-periodic. 

We check that (3(e) = 0(e) up to kif^x, for some k G Z. Therefore, we show that r(T(0), 0) = 2kn 
for some k G Z. 

By the definition of -D we get that D(t, 0)x® = x°(t) = x(t,xo)e~ Xot xo or 

D(t, 0) = e- Xot D(0, 0)e 9l ^ Q , (6.52) 
with gi sufficiently smooth such that g\(0) = 0. 



We show that gi(T(0)) = 2kir for some k G Z and r(T(0),0) = ffi(T(0)). 

Using (IQol) . it follows that D(t + T(0),0) = D(t, Q)e g ^ Q . Since 5(0, 0) = 0, we get D(T(0),0) = 
D(0, 0) = e- XoT( -^D(0, 0) e si( T ( ))O. Taking into account that |X(0)| T(0) = 2vr, we get ^ T ^))Q = 
I3. It follows that <7i(T(0)) = 2kn for some k G Z. 

Also, #(i,C£,0) = Z?(t,0)e- T (*'°) Q = e - Xoi Z)(0, O^WQe-^*' ^ by (lfi~52l . On the other hand, 
*(£,C \0) = e- x °*D(0,0). Then, we get r(i,0) = 51 (t) + 2/(t)vr with Z(i) G Z. Since 51 and r are 
sufficiently smooth, we get that l(t) is sufficiently smooth and, since l(t) G Z, we get l(t) =constant. 
Since r(0, 0) = 51(0) = 0, we get I = 0. 
Thus r(t,0) = (t) implies r(T(0),0) = ffi(T(0)) = 2/cvr. 

Then, /?(0) - - _ 2M 



T(0) ~~ T(0) ~~ T(0)- 



Therefore, (3(e) = Mjfc + O(e) = A;^ + /3i(e), where /?i(e) = O(e). 



We have C MRW (t,e) = B*(t,e)e~^ e ^ t = B*(t,e)e'^ kQt e -^ £ ^ t . 

Let B1(t,e) = B*(t,e)e~^ kQt , B^(t,e) is T(e)-periodic and BJ(0,e) = 5*(0,e). 
We may drop the subscript 1 in B\ and (3\(e). 

The stability issue for C MRW (t, e) is proved in 0. It is now easy to get the results for the differen- 
tial equations (|3.3j) using the fact that *(C, t) = [&(A, t)] , where A = C . This ends the proof 
of the third conclusion of Proposition 15.11 and its proof. □ 

Proof of Theorem \5. 6 A It is a consequence of Theorem 13.21 Proposition 15.11 and of the SO(2)- 
equivariance of the flow 

□ 
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